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Abstract 

We study theoretical implications of direct dark matter searches in the minimal super- 

Q^ symmetric standard model (MSSM). We assume that no accidental cancellations occur 

A in the spin-independent elastic neutralino-quark scattering cross section, but do not im- 

OJ pose any relations among the weak-scale MSSM parameters. We show that direct detec- 

r^, tion cross section below 10~^^ cm^ requires the lightest supersymmetric particle (LSP) 

neutralino to be close to either a pure gaugino or pure Higgsino limit, with smaller cross 

sections correlated with smaller admixture of the subdominant components. The cur- 

QQ rent XENONIOO bound rules out essentially all models in which the lightest neutralino 

"^ has the Higgsino fraction between 0.2 and 0.8. Furthermore, smaller direct detection 

cross sections correlate with stronger fine-tuning in the electroweak symmetry breaking 

sector. In the gaugino LSP scenario, the current XENONIOO bound already implies 

^ some fine-tuning: for example, at least 10% tuning is required if the LSP mass is above 

^-H 70 GeV. In both gaugino and Higgsino LSP scenarios, the direct dark matter searches 

currently being conducted and designed should lead to a discovery if no accidental 

^ cancellations or fine-tuning at a level below 1% is present. 
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1 Motivation and Philosophy 

The Minimal Supersymmetric Standard Model (MSSM) is a well-motivated and extensively 
studied extension of the Standard Model (SM) at the electroweak scale. Originally introduced 
as a way to solve the gauge hierarchy problem of the SM, the model turned out to have many 
other intriguing features. In particular, the MSSM contains an attractive particle dark matter 
candidate, the lightest neutralino [Tl[2]. If such a neutralino is the lightest supersymmetric 
particle (LSP) and R-parity is conserved, it is stable, and its thermal abundance naturally 
falls in the range indicated by the observed dark matter density. 

Probing the microscopic nature of dark matter requires observing interactions of individ- 
ual dark matter particles with ordinary matter via scattering or annihilation. An extensive 
experimental effort toward this goal is currently under way. In particular, direct detection 
techniques, which attempt to detect collisions between ambient dark matter particles and 
nuclei in the target, have seen significant advances in sensitivity with experiments such as 
CDMS 13], EDELWEISS [4J, and most recently XENONIOO [SIE]. So far, the results of these 
searches are null: no evidence for dark matter has been observedj^ The goal of this paper 
is to investigate the implications of these results, and improvements that are likely to come 
in the next few years, for the MSSM, with the assumption that the dark matter is made 
entirely of neutralinos. 

In the standard minimal WIMP framework, which is applicable to the MSSM dark matter 
throughout the model parameter space, the results of direct detection searches are presented 
as upper bounds on (or perhaps, in the future, measurements of) the spin-independent 
neutralino-proton elastic scattering cross section at zero momentum exchange Jj which we 
will call "direct detection cross section" for short. As usual with the MSSM, the theoretical 
prediction for this cross section depends on a large number of model parameters. The tradi- 
tional approach is to reduce the number of parameters by assuming relations among them, 
typically at the high energy scale {e.g. mSUGRA/cMSSM). Even then, the cross section 
predictions vary widely depending on the parameters, and the results are usually presented 
as scatter plots resulting from scanning the parameters within some broad ranges. While 
such plots provide a useful target for experiments, much potentially interesting information 
is missing. A typical scatter plot shows the cross section varying over several orders of mag- 
nitude, and it is not clear what features of the model correspond to points with higher, or 
lower, cross sections. Can any qualitative statements about the MSSM be made given the 
cross section bound of 10~^^ cm^ (roughly corresponding to the current XENONIOO result). 



^DAMA [7 and CoGeNT [8 collaborations observed effects inconsistent with known backgrounds, which 
may be due to dark matter scattering. If so, the MSSM dark matter is disfavored, since it falls under 
the minimal WIMP framework which does not provide a good simultaneous fit to these experiments and 
XENONIOO [6 . At the moment, however, the experimental situation is quite confusing, and we will not 
take DAMA and CoGeNT into account in this study. 

^Converting the actually measured event rates into a cross-section bound or measurement requires a num- 
ber of assumptions, such as local dark matter density (see e.g.p^) and velocity distributions (see e.g.[TT), 
isospin symmetry of WIMP-nucleon couplings (see e.g.[15 ), etc. We will not discuss the potential uncer- 
tainties introduced by the conversion in this paper. 
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Figure 1: Feynman diagrams contributing to spin-independent elastic scattering of the neu- 
tralino dark matter particle on a nucleon in the MSSM. 



or some lower future bound? This will be the main focus of this paper. In particular, we 
will demonstrate a correlation between the direct detection cross section and the amount 
of fine-tuning in the electroweak sector: roughly speaking, model points with lower direct 
detection cross sections are more fine-tuned. 

In order to be as general as possible, we will treat all weak-scale MSSM parameters as in- 
dependent, without assuming any relations among them. The tree-level processes contribut- 
ing to the direct detection cross section are shown in Fig. [T] The key assumption underlying 
our analysis is that no accidental cancellations take place among various contributions to 
the direct detection cross section in the MSSM. By "accidental", we mean a cancellation 
which is exact only on a measure-zero hypersurface inside the full MSSM parameter space. 
Equivalently, an accidental cancellation is indicated by an anomalous sensitivity of the cross 
section to the MSSM parameters (measured, for example, by its logarithmic derivative) along 
at least one direction in the parameter space. In particular, any cancellation between the s- 
and t-channel diagrams in Fig. [T] would be accidental, since they depend on diflFerent sets of 
MSSM parametersjj Thus, for making qualitative statements, it is sufficient to consider only 
one of the diagram classes; the other one will, at worst, produce an order-one correction to 
the cross section. We will focus on the t-channel diagrams. Fig. [l] (a). We make this choice 
because three of the five MSSM parameters which enter these diagrams, /i, tan/3, and m^, 
also enter the tree-level prediction for the Z mass. In this way, the direct detection cross 
section is connected to electroweak symmetry breaking. 

When comparing with experimental data, we will assume that the local dark matter 



^Of course, different MSSM parameters may be related once the SUSY-breaking sector is understood, 
so that a cancellation that appears accidental from the weak-scale point of view may in fact be natural in 
the full theory. Such a situation, however, appears extremely unlikely in the particular situations where we 
apply the "no accidental cancellation rule" in this study. For example, a cancellation between the s- and 
t-channel diagrams in Fig. [T] would require a complicated relation involving squark and gaugino soft masses, 
the /i parameter, tan/3, and the Higgs mass terms. It is very difficult to imagine a SUSY-breaking model 
producing such a relation. 



density is at its canonical value, 0.3 GeV/cm^, and that dark matter is made entirely of the 
MSSM neutralinos. However, the main results of the paper concern theoretical predictions 
for direct detection cross section and are independent of this assumption. Moreover, in 
most of the analysis, we will not impose the constraint that the neutralino relic density 
predicted by the standard thermal decoupling calculation be consistent with observations. 
This is motivated partly by the desire to keep the analysis as general as possible: Statements 
made without the relic density constraint would be applicable even if there is non-thermal 
production of neutralinos, late inflation, other deviations from standard FRW cosmological 
evolution, etc. More technically, it is motivated by the fact that the object we focus on, 
the t-channel direct detection cross section, depends on just flve MSSM parameters, while 
the relic density is a function of many more. Imposing this constraint would thus greatly 
complicate the analysis. (As an exception, we will impose a mild version of the relic density 
constraint in the part of the analysis dealing with Higgsino dark matter, where this will be 
required to make an interesting statement about flne-tuning.) 

Before we proceed, let us remark on a few related analyses in the literature. The con- 
nection between direct detection cross section and flne-tuning in the electroweak sector was 
noted, in the mSUGRA context, in Ref. j9] (see also [ID]). A study of this connection in 
a general MSSM framework similar to ours, focusing on the bino-LSP scenario, appeared 
in Ref. [llj. A recent study of "generic" direct detection cross sections in a general MSSM 
framework was presented in Ref. [I2]. While our approach is similar, our setup is even more 
general: in particular, we do not require gaugino mass uniflcation, and (for the most part) do 
not impose the thermal relic density constraint. Electroweak flne-tuning, which did not enter 
the analysis of [12j, plays the central role in our discussion. We also discuss the impact of 
the recent XENONIOO bound in the general MSSM framework, with the main assumption 
being the absence of accidental cancellations. This discussion complements Refs. fT6l fTTj 
which interpreted the XENONIOO result within mSUGRA and other constrained MSSM 
frameworks. 

The rest of the paper is organized as follows. In section [2} we set the notation and collect 
the main formulas used in the analysis. The main results of the paper are contained in sec- 
tion |3| where we present a set of scatter plots demonstrating correlations between the direct 
detection cross section and physical quantities such as Higgsino fraction of the neutralino 
and the amount of flne-tuning in the electroweak symmetry breaking. These plots assume 
real, positive MSSM parameters. This assumption is lifted in section |4| where negative and 
complex soft masses are considered, and it is shown that the interesting correlations persist 
once the points with accidental cancellations in the cross section are eliminated. Finally, we 
conclude in section [H 



2 Analysis Setup 

The direct detection cross section has the form [2] 
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where m^ is the neutrahno-proton reduced mass and 



. = M5, (1) 



^ = E f^>, + I /S E A,. (2) 

P q=u,d,s q=c,b,t 

The nuclear formfactors are defined by 

f^: = ^^=^ , (3) 

TTZp 

and 

A'i = i- E /^^- (4) 

q=u^d^s 



In our numerical work, we will use the following values 

#2=0.08, /S= 0.037, /g=0.34. (5) 

It is well known that there is a significant uncertainty on these formfactors, in particular 
fs [18j. The uncertainty could easily be incorporated in our study; however, since our main 
interest is in qualitative trends rather than precise quantitative statements, we will not 
do so here. The dependence of the cross section on the underlying particle physics model 
is contained in the coefficients Aq^ which in the MSSM at tree level are computed with the 
Feynman diagrams in Fig.[T} As explained above, we will focus on the t-channel contribution, 
given by [T9] 
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+C^D^ (^ - -yi Re [5u{9Z^2 - g'Z^^)]] , (6) 

where for up-type quarks 

Bu = sin /3 , Cu = sin a, Du = cos a, 5iu = Z^s , 52u = Z^4 ; (7) 

while for down-type quarks 

Bd = cos /5 , Cd = C08 a, Dd = - sin a, 5id = Z^4 , d2d = -Z^s . (8) 



The notation here is consistent with that used, for example, in Martin's review [20]; in 
particular, the lightest neutralino is 

XS = Z^rB + Z^,W' + Z^sH', + Z^4^° , (9) 

and h and H are the two CP-even Higgs bosons {rrih < ttih)- The coefficients Ai can be 
expressed in terms of only five MSSM parameters: 

p, = (Mi,M2,/x,tan;5,mA), (10) 

where Mi, M2, and jjl appear in Eq. ^ implicitly via the neutralino mixing matrix elements. 
These parameters will always be defined at the weak scale; since no unification or any other 
relation among the parameters is assumed, we do not need to consider their renormalization 
group evolution. In general, the parameters pi are complex; however it can be shown (see, 
for example, Ref. [30j) that only two phases are physical: 

Lpi = arg(/iMi sin2^), Lp2 = arg(/xM2sin2;5) . (11) 

These phases are constrained by measurements of electric dipole moments (EDMs) [21], 
although maximal phases are allowed if squarks and sleptons are very heavy [22j . The light 



Higgs mass rrih can be expressed in terms of the parameters in (10) at tree level. It is of 
course well known that a large loop correction is required to satisfy the LEP-2 lower bound 
on rrih] this correction is dominated by the top and stop loops and including it would bring 
in a few additional MSSM parameters into the game. In this study, we avoid doing this by 
simply fixing rrih at a fixed value consistent with LEP-2, rrih = 120 GeV. (In the MSSM, the 
upper bound on rrih is about 135 GeV; a variation of ruh in the allowed range does not have 
a strong effect on the direct detection cross section.) In other words, we assume that for any 
set of p^'s, the other MSSM parameters can be chosen so that rrih = 120 GeV. 

The other quantity of interest for us is the amount of fine-tuning in the electroweak 
symmetry breaking (EWSB) sector of the model. A tree-level analysis will suffice. (The 
discussion below is taken from Ref. [23].) The key formula is the relation of the Z boson 
mass to the MSSM Lagrangian parameters: 



ml = -mill ]-nil{l + ^—]-2\^i\\ (12) 




where m^ and m^ are the Lagrangian masses for the up- type and down- type Higgs doublets, 
and 

sm2/3= ^^ (13) 

We quantify fine-tuning by computing 
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6, ji are the relevant Lagrangian parameters. (This is analogous to the 



fine-tuning measure introduced by Barbieri and Guidice [24j, although here it is applied to 
weak-scale, rather than Planck/GUT-scale, MSSM parameters.) Using the well-known tree- 
level relations to express m\ and m^ in terms of the parameters listed in (10), we obtain 
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where we assumed tan/5 > 1. The overall fine-tuning A is defined by adding the four 8's 
in quadruture; values of A far above 1 indicate fine-tuning. The qualitative behavior of 
fine-tuning is easy to understand by taking the limit of large tan/5, where the parameters 
5{m\) and 5(m^) are small, and 
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Thus, increasing /x, and to a lesser extent m^, requires fine-tuning. On the other hand, as 
(3 approaches 7r/4, the factors of 1/ cos 2^5 and tan 2^, present in all four 5's, become large, 
and as a result the model is always fine-tuned for tan/5 ^ 2. 

Our strategy for the rest of the paper is as follows: We will perform scans over the five 
parameters p^, computing the direct detection cross section and the fine-tuning measure 
discussed above for each point, and study the correlation between these two quantities as 
well as with other relevant parameters such as the LSP mass, Higgsino/gaugino fractions, 
etc. We will then explain physical reasons for each observed interesting correlation. 

3 Results: Real, Positive Parameters 

We begin by performing a restricted scan in which we assume that all five pi parameters 
are real and positive. All the correlations that we found show up most clearly in this scan, 
making it a good place to start the discussion. In more complete scans, the correlations 
persist, but are somewhat obscured by the possibility of accidental cancellations within the 
i-channel contribution to the direct detection cross section. This will be discussed in detail 
in the next section. 
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Figure 2: Left panel: Direct detection cross section vs. Higgsino fraction of the neu- 
tralino. Right paneh Direct detection cross section vs. the dark matter particle mass, 
for points with purity above 0.2 (red), between 0.1 and 0.2 (orange), 0.01 and 0.1 (green), 
and 10~^ and 0.01 (cyan). The lines correspond to the XENONIOO 100 days exclusion 
limit [6] (black/solid), and the projected sensitivities of the XENONIOO upgrade [25l [26] 
(blue/dotted) and XENON-IT ^26j (red/dashed). Real, positive values of the scanned MSSM 
parameters are assumed. 



A set of 10^ MSSM points was generated, distributed randomly, uniformly in log Mi, 
logM2, log/x, logm^, and tan^S, within the following scan boundaries: 



Ml G [10, 10^] GeV; M2 G [80, 10^] GeV; 
/i G [80, 10^] GeV; tua G [100, 10^] GeV; 
tan;5G [2,50]. 



(17) 



We compute the neutralino and chargino masses for each point in the scan, and exclude 
points with at least one chargino below 100 GeV (excluded by LEP-2), as well as those 
where the lightest chargino mass is below the lightest neutralino mass. We do not impose 
any other experimental constraints, since they depend on the MSSM parameters beyond the 
five Pi being scanned here, and thus generically can be satisfied by varying those parameters 
for any given pi. The scatter plots in this section are based on the 73064 points that pass 
these constraints. 

3.1 Higgsino Fraction Constraint 

Our first result concerns the correlation between the direct detection cross section and the 
Higgsino fraction of the neutralino, defined as 



F, 
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\Z^s\ + |Z 



X4| 



(18) 



As is clear in Fig. [2| a direct detection cross section limit below (a few)xlO~^^ cm^ puts a 
constraint on the Higgsino fraction, requiring that it be close to either (the "pure gaugino" 
case) or 1 (the "pure Higgsino" case). It is convenient to define "neutralino purity" p as 

j9 = min(F^,l-F^). (19) 

The bound placed by the XENONIOO experiment [6j already rules out essentially all MSSM 
dark matter models with p > 0.2, and most models with p > 0.1, especially for the LSP mass 
above 50 GeV.QJThe proposed XENONIOO upgrade [251 [26] will be able to probe values of p 
down to 0.01 for the LSP masses above 50 GeV, while a 1-ton upgrade [26] will have a reach 
down to p ^ 10~^ through most of the mass rangejj If the dark matter is not discovered 
at that stage, the only possibility in the MSSM would be a pure gaugino or Higgsino with 
< 0.1% admixture of the other components. 

The physical origin of this constraint is easy to understand. In the gauge eigenbasis for 
neutralinos, the neutralino-neutralino-Higgs couplings have the form [gW^ -\- g' B)H H ] there 
are no gaugino-gaugino-Higgs or Higgsino-Higgsino-Higgs couplings in the MSSM. In the 
mass basis, the couplings have the form 

X^X^h : {gZ^2 - g' Z^i){cosaZ^^ + sinaZ^g) , 
X^X^H : {gZ^2 - g' Z^i){smaZ^^ - cosaZ^a) . (20) 

If the x^x^h is of its natural size {i.e. no accidental cancellations or small mixing angles 
are present), the direct detection cross section from t-channel Higgs exchange is of order (a 
few)xlO~^'^ cm^ or above. Barring accidental cancellations, the only way to obtain a smaller 
cross section is to suppress this coupling by choosing the LSP to be an almost pure gaugino 
or Higgsino, which is precisely what is seen in Fig. [2j 

3.2 Gaugino Dark Matter and Elect roweak Fine- Tuning 

To discuss the connection between direct detection cross section and EWSB fine-tuning, it is 
useful to divide the scan points into two sets: the points where ji < Mi and fi < M2, and the 
rest. We will refer to the first set of points as the "Higgsino LSP" sample, while the second 
set will be called "gaugino LSP" sample. Of course, these names are not precise, since each 
sample contains points with /x ^ Mi^2 where the LSP is a roughly equal mixture of the two; 
however, as we saw above, such points always have high direct detection cross sections and 
will not infiuence the fine-tuning discussion. (Most of these points are in any case already 
ruled out by XENONIOO, though for simplicity we will not impose this constraint here.) 



^A well-known example of such a model is the "well-tempered neutralino" scenario [30 . The fact that 
this scenario is disfavored by XENONIOO has already been noted in Ref. [16j . 

^Many next-generation direct dark matter searches have been proposed, such as XMASS [27 , LUX [28] . 
and superCDMS [29 . Needless to say, we use projections from the XENON collaboration simply as a 
benchmark, and do not mean to endorse or express a preference for a particular technology or experimental 
proposal. Projected sensitivities of any proposed experiment can be easily superimposed on our plots. 
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Figure 3: Direct detection cross section vs. fine-tuning of electroweak symmetry breaking, for 
gaugino-like neutralino (real, positive values of p^, i = 1 . . . 5, are assumed). Cyan, green and 
red points correspond to the dark matter particle mass in the intervals [10, 100), [100, 1000), 
and [1000, 10^] GeV, respectively. The cyan, green and red lines show the analytic lower 
bound in Eq. (21), with Mlsp = 10, 100, 1000 GeV, respectively. Real, positive values of the 



MSSM parameters are assumed. 



Let us first consider the gaugino LSP sample. The correlation between the direct detec- 
tion cross section and the amount of fine-tuning in the EWSB sector (quantified by A defined 
in Sec. [2]) in this sample is shown in Fig. [3] For a given amount of fine-tuning, the direct 
detection cross section cannot be reduced below a certain lower bound, with smaller cross 
sections possible only for more finely-tuned models. The physical origin of this correlation 
is again simple: as we saw above, the cross section can only be reduced by reducing the 
Higgsino admixture in the LSP, and the only way to achieve this is to raise fi. But doing 



so increases the fine-tuning, as is easily seen from Eq. (12) or Eq. (16). It should also be 



mentioned that the sparsely populated region in the top right corner is merely a consequence 
of our choice of scan boundaries and distribution (points in this region require both large /x 
to be greatly fine-tuned and comparably large Mi or M2 to have large cross sections, which 
is extremely unlikely given our scan points are log distributed in these parameters). 

Also plotted in Fig. |3] is an analytic expression for the minimal direct detection cross 
cross section possible for a given amount of fine-tuning. For a given tan/3 and LSP mass, it 



is given by (for derivation, see Appendix A): 
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It is clear that the lowest possible cross section for fixed A requires the highest possible tan (3 
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Figure 4: Left panel: Direct detection cross section vs. dark matter particle mass, for 
gaugino-like neutralino. Red, green and cyan points correspond to EWSB fine-tuning in the 
intervals (0,10); [10,100); [100,1000), respectively. Right panel: Minimal direct detection 
cross section compatible with EWSB fine-tuning of 10 (red line), 100 (green line), and 1000 



(cyan line) according to Eq. (21). Current and projected XENON bounds are superimposed 
(notation as in Fig. [2]). Real, positive values of the MSSM parameters are assumed. 



and the lowest possible Mlsp- In the plot in Fig. [3| we used tan^S = 50, corresponding to the 
upper boundary of the scan. Another noteworthy feature is that the lowest possible direct 
detection cross sections occur in the Higgs decoupling limit, vrtA ^ ^z- 

The practical implication of the correlation in Fig.[3]is clear: in the gaugino LSP scenario, 
a sufficiently strong bound on direct detection cross section implies a non-trivial fine-tuning 
in the EWSB sector of the MSSM. The precise level where the cross section bounds become 
relevant for fine-tuning depends on the LSP mass. This is illustrated in Fig. |4} The current 
XENONIOO bound is already relevant: for example, fine-tuning of 1 part in 10 or less is 
only possible for LSP masses below about 70 GeV. Null result of the proposed XENONIOO 
upgrade [25l [26] would imply at least 1% fine-tuning for WIMPs at 100 GeV or above. 
This level of fine-tuning is similar to what is required to accommodate the non-observation 
of the Higgs boson at LEP-2, the famous "little hierarchy problem" j3l] of the MSSM. A 
XENON-IT upgrade [26j would be able to probe all models fine-tuned at 1% level or better, 
and most models fine-tuned at 0.1% or better. No observation at the level of sensitivity 
where irreducible neutrino backgrounds would render further progress impossible, about 
10~^^ cm^ [32j, would imply fine-tuning of at least 0.1% across the entire LSP mass range. 
We stress that these statements are remarkably general: we only assumed that a mostly- 
gaugino LSP contributes all of the present dark matter, and that no accidental cancellations 
occur in the direct detection cross section. No assumptions whatsoever have been made 
about the SUSY-breaking model, and the neutralino thermal relic density constraint has not 
been imposed, so that the result is independent of cosmological history. (The discussion here 



10 



also assumes that p^'s are real and positive, but this restriction does not substantially affect 
the results, as shown in the next section.) 

3.3 Higgsino Dark Matter 

Now, let us switch our attention to the Higgsino LSP sample. It is clear that, if no additional 
assumptions are made, no correlation between direct detection cross section and fine-tuning 
exists in this case. Indeed, if /x ^ 100 GeV, no fine-tuning is necessary in the EWSB, while 
the direct detection cross section can be suppressed by choosing Mi^2 ^ M- Large values of 
Mi^2 only affect fine-tuning at the one-loop level, and these parameters can be at the multi- 
TeV scale without significant tuning. In this situation, the direct detection cross sections are 
many orders of magnitude below the current sensitivity, and can even be below the 10~^^ 
cm^ level where the neutrino background would render direct detection impossible. 

To make an interesting statement in this situation, we need to make an additional as- 
sumption. An obvious one is to demand that the LSP has the correct relic density, assuming 
conventional FRW cosmology and no non-thermal production. The relic density is typically 
determined by the LSP annihilation cross section. In the limit of pure Higgsino LSP, the 
next-to-lightest neutralino and the lightest chargino are quasi-degenerate with the LSP, and 
co-annihilations among these states need to be taken into account in the relic density cal- 
culation [34]. The cross sections are typically dominated by annihilations into electroweak 
gauge bosons, Vl^'s and Z's. At tree level, these cross sections are completely determined 



by the same five MSSM parameters, Eq. (10), that entered our analysis of direct detection. 



There are, of course, other contributions to the annihilation cross section, such as the dia- 
grams with top final states. Including these processes would introduce more parameters and 
complexity into our analysis. To avoid this, while still keeping open the possibility that they 
contribute significantly to the cross section, we impose a one-sided relic density constraint: 
we demand that the cross section of (co) annihilation into W/Z final states be no larger thdiH 
what is required to obtain the observed relic density, ^dmh^ = 0.110 ± 0.006 [33]. 

For each point in the Higgsino LSP sample, we performed the relic density calculation 
using the numerical package DarkSUSY [35j. In DarkSUSY runs, all mass scales other 
than those in the five input parameters in (fTol) were set to 10 TeV, effectively eliminating 
annihilations into quarks and leptons. We require that the relic density calculated with 
DarkSUSY not be lower than 2a below the observed relic density. Once this constraint is 
imposed, a simple relation between fine-tuning and direct detection cross section emerges, 
seen clearly in Fig. [5] If direct detection cross section is constrained to be below about 



2 X 10 ^^ cm^, the LSP must be a pure Higgsino, as we saw in section 3.L For pure 



Higgsino, the annihilation cross section is too large, unless /x ^ 1 TeV. But such large values 



of /i require fine-tuing in the EWSB of at least about 0.25%, as is easily seen from Eq. (16). 
Thus, a direct detection bound of about 10""^"^ cm^ for a 1 TeV LSP mass would imply a 
"little hierarchy problem" for Higgsino LSP. The current XENONIOO bound at 1 TeV mass 
is about 8 x 10""^^ cm^, so no such statement can yet be made. The proposed XENONIOO 
upgrade [25l[26] should reach the required sensitivity. If no signal is discovered, the Higgsino 
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Figure 5: Direct detection cross section vs. fine-tuning of electroweak symmetry breaking, 
for higgsino-like neutralino. Real, positive values of the MSSM parameters are assumed, and 
a one-sided thermal relic density constraint, ^pred ^ ^obs, is imposed, as explained in the 
text. 



LSP scenario would be inconsistent with naturalness of EWSB at a sub-percent level. 

4 Results: Full Scan 

In this section, we remove the requirement of real, positive pi. Since complex phases of the 
MSSM parameters are generally constrained by measurements of EDMs, we first remove the 



condition pi > 0, but keep them real. Only two phases are physical, see Eq. (11), and we 
choose the basis where /x and the Higgs vevs are positive but Mi and M2 can have either 
sign. We generate a set of 10^ MSSM points distributed randomly, uniformly in log|Mi|, 
log IM2I, log/i, logm^, and tan/3, within the following scan boundaries: 



|Mi| G [10, 10^] GeV; IM2I G [80, 10^] GeV; 
/i G [80, 10^] GeV; tua G [100, 10^] GeV; 
tan/3G [2,50]. 



(22) 



The signs of Mi and M2 are chosen between + and — with equal probability. After imposing 
the same requirements as in the positive-only scan of section [3| we are left with 76546 points, 
which are included in the scatter plots below. 

Fig. |6] demonstrates the correlation between the direct detection cross section and the 
Higgsino fraction of the neutralino in this sample. The bulk of points show a correlation 
similar to that in the real-only sample: Most points with very small direct-detection cross 
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Figure 6: Direct detection cross section vs. Higgsino fraction of the neutralino. Green, orange 
and red points correspond to Aacc below 10, between 10 and 30, and above 30, respectively. 



sections are close to either pure-gaugino or pure-Higgsino limit. However, there are some 
"outliers" , which have direct detection cross sections well below 10~^^ cm^ in spite of having 
order-one Higgsino and gaugino fractions. The reason for this is accidental cancellations 
within the t-channel contribution to the cross section. An accidental cancellation in the 
cross section at a particular point in the parameter space is characterized by its anomalous 
sensitivity to the input Lagrangian parameters at that point. To quantify this, we introduce 
the measure 



\ 



E 



' d log a 
d log Pi 



(23) 



where pi = (Mi, M2,/x, tan/5,mA). The points in Fig. |6| are color-coded according to this 
measure, making it clear that the outlier points are uniformly characterized by severe ac- 
cidental cancellations. Thus, the conclusion of the analysis in the previous section remains 
unchanged: Limits below 10""^^ cm^ imply that the MSSM dark matter neutralino has to 
be either pure gaugino or pure Higgsino, unless there are accidental cancellations in the 
direct detection cross section. This is further illustrated by Fig. [7| No firm conclusions on 
the Higgsino fraction can be drawn from the current or projected XENONIOO bounds if 
all points in the scan are included (left panel). However, once the points with significant 
accidental cancellations are excluded (right panel), the situation becomes similar to the case 
of positive-only parameters (compare with Fig. [2]), and generic order-one Higgsino-gaugino 
mixture is already strongly disfavored. 

A similar picture emerges for the correlation between direct detection cross section and 



EWSB fine-tuning in the gaugino LSP case (see Sec. 3.2). Fig. |8| shows this correlation. The 
vast majority of points in the scan obey the lower bound on the cross section as a function 
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Figure 7: Left panel: Direct detection cross section vs. the dark matter particle mass, for 
points with purity above 0.2 (red), between 0.1 and 0.2 (orange), 0.01 and 0.1 (green), and 
10~^ and 0.01 (cyan). Current and projected XENON bounds are superimposed (notation 
as in Fig. [21). Right panel: Same, including only points with Aacc ^ 5. 



of A, Eq. (21). Most points that violate the bound, and all that violate it by an order of 



magnitude or more, suffer from severe accidental cancellations. In fact, a slightly stronger 
version of the analytic lower bound. 



:i.2 X 10"^^ cm^) 



120GeV\^ 1 



A 



rrih 



mm 



1 1 Mlsp 



tanks' v^ mz 



(24) 



is obeyed by all but a few of models with no strong accidental cancellations. 

The direct detection/EWSB fine-tuning correlation is further illustrated by Fig. |9J If all 
points are included (left panel), the correlation is somewhat washed out, with a fairly large 
number of "natural" MSSM points yielding cross sections well below current or projected 
XENON bounds, and even some points with cross sections below the neutrino background of 
10""^^ cm^ [i32j. Most of these points, however, are characterized by severe accidental cancel- 
lations, and once they are removed (right panel), the correlation between cross section and 
EWSB tuning is almost as clear as in the positive-only case (compare with Fig. ^. Thus, 
once points with accidental cancellations are discarded, the discussion of fine-tuning implica- 



tions of cross section bounds at the end of Sec. |3.2| qualitatively applies to the more general 
case where MSSM soft masses of either sign are allowed, although quantitative statements 
are slightly weaker: for example, a few points with EWSB fine-tuning of 1 part in 10 or less 
and LSP masses in the 70-100 GeV range survive the current XENONIOO bound. 

A further demonstration of the accidental nature of cancellations leading to low direct- 



detection cross sections is given in Fig. [T0| where we fix three of the five MSSM parameters 
and scan over the other two. Points with low cross sections (orange, below 10""^^ cm^, and 
red, below 10~^^ cm^) are clearly seen to be confined to narrow bands within the parameter 
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Figure 8: Direct detection cross section vs. EWSB fine-tuning. Green, orange and red points 
correspond to Aacc below 10, between 10 and 30, and above 30, respectively. Also shown is 



the analytic lower bound in Eq. ( [2l| ) (blue/solid line), and its modified version in Eq. (24) 
(black/dashed line). 



space. 

Finally, we repeated the scan allowing Mi and M2 to have random phases between and 
2n. This may be allowed if squarks and slept ons are very heavy, suppressing contributions 
to EDMs ^^. The results of the scan are broadly similar to those for real parameters. 
Again, the correlations observed in section [3] persist, but are somewhat obscured by points 
with accidental cancellations. As an example, in Fig. [TT] we show a scatter plot of direct 
detection cross section vs. EWSB fine-tuning for models with a gaugino-like LSP. We use 
Aacc to quantify accidental cancellations, extending its definition to include the logarithmic 
derivatives with respect to the two phases in the sum. Again, the vast majority of points 
obey the lower bound on the cross section as a function of fine-tuning, and the ones that do 
not typically have Aacc ^ 1, indicating accidental cancellations. 



5 Conclusions 

In this paper, we considered the theoretical implications of direct dark matter searches in 
the MSSM framework. Unlike almost all previous analyses, we did not require any relations 
among the MSSM parameters, such as mSUGRA, and did not impose thermal relic den- 
sity constraints. Instead, we make the MSSM parameter space tractable by assuming that 
no accidental cancellations take place among the 5-channel and t-channel contributions to 
the spin-independent elastic neutralino-quark cross section, and focusing on the t-channel 
diagrams. These only depend on 5 MSSM parameters. We performed extensive scans over 
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Figure 9: Left panel: Direct detection cross section vs. dark matter particle mass, for 
gaugino-like neutralino. Red, green and cyan points correspond to EWSB fine-tuning in the 
intervals (0,10); [10,100); and[100, 1000], respectively. Right panel: Same, including only 
points with Aacc ^ 5. Current and projected XENON bounds are superimposed (notation 
as in Fig. [2]). 



these parameters, in order to identify correlations between the direct detection cross sections 
and other quantities of physical importance. The following simple picture emerges from our 
analysis: 



• If the LSP is a generic mixture of Higgsino and gauginos with order-one admixture 
of each component, the direct detection cross section is always above 2 x 10~^^ cm^. 
Such models are already severely constrained: for example, essentially all models with 
a Higgsino fraction between 0.2 and 0.8 are ruled out by XENONIOO. 

• Lowering direct detection cross section below the 10~^^ cm^ level requires that the 
LSP be either pure gaugino (bino or wino), or pure Higgsino. In both cases, smaller 
cross sections correlate with higher purity (i.e., smaller admixture of the subdominant 
components) of the LSP. 

• If the LSP is a gaugino, smaller direct detection cross sections correlate with stronger 
fine-tuning in the elect roweak symmetry breaking sector, since they require higher 
values of the /x parameter. The current XENONIOO bound already implies non-trivial 
fine-tuning, albeit rather mild at this point: for example, 1 part in 10 or worse tuning 
is required if the LSP mass is above 70 GeV. Future experiments could put stronger 
stress on the model: for example, if no signal is seen at the proposed XENONIOO 
upgrade, fine-tuning of 1% or worse would be required for LSP mass above 100 GeV. 
This would present a new "little hierarchy problem" for the MSSM (at least if the idea 
that the neutralino makes up all of the dark matter is taken seriously). 
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Figure 10: Left panel: Scatter plot of direct detection cross section as a function of tan^S 
and — M2, with Mi = —150 GeV, /x = 200 GeV, and ttia = 500 GeV. Cyan, green, orange 
and red points have log^g crcm2 > — 45, log^g crcm2 ^ (— 46, — 45),logiQacm2 ^ (—46,-47) and 
logiocrcm2 < —47, respectively. Right panel: Same, as a function of fi and —Mi, with 
M2 = -200 GeV, tan/3 = 10, and tua = 500 GeV. 



• If the LSP is a Higgsino, no correlation between direct detection cross section and 
fine-tuning can be established. However, if an additional mild assumption, the one- 



sided thermal relic density constraint (see sec. 3.3) is imposed, all points with direct 
detection cross section below about 2 x 10~^^ cm^ have an LSP mass of 1 TeV or above, 
and electroweak fine-tuning at the level of 0.25% or worse. Non-observation of a signal 
at the proposed XENON 100 upgrade would imply this level of tuning in the Higgsino 
LSP scenario. 

Once again, it is worth emphasizing that these conclusions are very general, and apply in 
the full "phenomenological" MSSM, without undue theoretical prejudice [36]. The degree of 
fine-tuning in the EWSB is a widely accepted quantitative "figure of merit" used to assess 
relative theoretical attractiveness of various regions of the MSSM parameter space. We 
established a clear correlation between this measure and the direct detection cross section. 
The main conclusion of our analysis is that, if the MSSM is the true model of microscopic 
physics and dark matter, and no fine-tuning at a sub-percent level is present, the direct dark 
matter searches currently being conducted and designed should lead to a discovery. 

A similar analysis could be performed in other models of electroweak symmetry breaking 
with particle dark matter candidates. An interesting direction for future work is the so-called 
next-to-minimal supersymmetric standard model (NMSSM), an MSSM with an additional 
singlet field in the Higgs sector. This model significantly alleviates the fine-tuning due to 
the Higgs mass lower bound, and the additional "singlino" admixture may be present in 
the LSP, affecting the direct detection cross section. It would be interesting to see if the 
correlations discussed here persist in the NMSSM. 



17 



10-44 
10-45 



10-48 
10 



49 






^ 







10 



50 100 500 1000 

fine tuning 



5000 



Figure 11: Direct detection cross section vs. EWSB fine-tuning, for the scan with arbitrary 
phases of Mi and M2. The notation is identical to Fig. [8| 



Acknowledgments 

This research is supported by the U.S. National Science Foundation through grant PHY- 
0757868 and CAREER grant No. PHY-0844667. This research was supported in part by 
the National Science Foundation under Grant No. PHY05-51164. MP is grateful to Michael 
Peskin and to Aaron Pierce for useful remarks. MP acknowledges the hospitality of the 
Kavli Institute for Theoretical Physics at Santa Barbara, where part of this work has been 
completed. BS acknowledges the hospitality of the Theoretical Advanced Study Institute 
(TASI-11) at the Univeristy of Colorado at Boulder. 



A Analytic Direct Detection/Fine- Tuning Relation for 
Gaugino LSP 



In this appendix, we will derive the formula (21) for the minimal direct detection cross 



section consistent with fixed values of tan/3, Mlsp, and A. We assume real and positive /x. 
Ml and M2, as in the partial scan of section [3| and focus on the gaugino LSP region, defined 
by /x > Ml and/or /x > M2. As discussed in section [3| a small direct detection cross section 
requires a pure gaugino LSP, which implies /x ^ Mi and/or // S> M2, and Z^3,Z^4 <C 1. 



Moreover, it is clear from Eq. (20) that the direct detection cross section is minimized if the 



LSP is predominantly bino, since g' < g. Motivated by these considerations, consider the 
limit 

Ml < /x < M2 . (25) 
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M^o 



(26) 



Ignoring terms suppressed by M2, the neutralino mixing angles are obtained by diagonalizing 

the mass matrix 

Ml -Cfss^mz Sfss^mz \ 
-Cf3S^mz -/i 

S/sSwmz -/i 

where sq = sinO^ cq = cos0, and subscript w denotes the Weinberg angle. Assuming rriz <C // 
(which is always the case in the region of interest) and working to second order in l//x, we 
obtain 



J 



^X3 



^X4 



mz f Ml 

Sw Sf3 + Cf3 



mz 



-C/3 - Sp- 






(27) 



Barring accidental cancellations, the direct-detection cross section is minimized in the Higgs 
decoupling limit, rriA ^ trz- In this limit, the well-known relationships 



sin 2a 
sin 2/3 



m\ + m\ 



-^ -1 



tan 2a m^ 



m% 



m 



H 



mt 



tan 2/3 



-^ 1 



m. 



rrir^ 



(28) 



imply sino; ~ —cos/?, cos a ~ sin/5. Plugging this and Eq. (27) into Eqs. (J6]), ([T]), and ^ 
and ignoring the heavy Higgs exchange diagrams, yieldqj 



A, 



Aa 



wa 1 



sin 2/3 + 



Mlsp' 



(29) 



where a is the fine structure constant, and we used Mlsp = Mi + 0(l//x). The final step 
is to relate jjl to the fine-tuning measure A. For jjl > Q^ the two terms in the brackets in 



Eq. (29) cannot cancel, and the direct detection cross section is minimized (for fixed /x and 
^Lsp) at large tan/3, since sin 2/3 ^ 2tan~-^/3 ^ in this limit. As discussed in section [2| 
the fine-tuning in this limit is dominated by 5(/x) and 5(6), see Eq. (16), so that 



rur 



'\ 



/|4 + 



m. 



tan^/5 



(30) 



It is clear that the smallest direct detection cross section for a given value of A is achieved 
for the largest possible jjl consistent with that A, which occurs when ttia <C Vtan^/x. (It can 
be easily checked that there are always values of ttia consistent with this condition which 
are still large enough to ignore the H exchange contributions to direct detection matrix 
elements.) In this situation. 






(31) 



^The approximate expression for the direct detection cross section in the bino-LSP limit also appeared 
inRef. [n]. 
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so that we can rewrite Eq. (29) in its final form 

47ra 1/1 ,1 Mlsp\ 



Au-A,^ 2 ^x\^-a + ^x-^ • (32) 

"^mzml VA ytsiup ^/A mz J 



Squaring and multiplying by the appropriate nuclear formfactors yields Eq. (21). 
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